Local Density of States for Individual Energy Levels in Finite Quantum Wires 
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The local density of states (LDOS) in finite quantum wires is calculated as a function of dis- 
crete energies and position along the wire. By using a combination of numerical density matrix 
renormalization group (DMRG) calculations and analytical bosonization techniques it is possible to 
obtain a good understanding of the local spectral weights along the wire in terms of the underlying 
many-body excitations. 
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There has been increasing interest in quasi one- 
dimensional quantum wires over the last few decades. 
Tunneling experiments have revealed signatures of strong 
correlations in different setups ranging from a powerlaw 
behavior of the temperature dependent conductance in 
carbon nanotubes [l], Q to being able to map out the 
momentum resolved tunneling in parallel quantum wires 
in a GaAs/AlGaAs heterostructure [3, 4, 5]. By now the 
local tunneling density near nanotube edges can be mea- 
sured fully resolved in energy and position [6] , which has 
been interpreted as the interference pattern of collective 
excitations. Standing waves at discrete energies corre- 
sponding to the few lowest lying levels have also been 
observed in finite tubes with screened interactions [3, H[ • 
It can be expected that similar experiments will be able 
to identify collective excitations in short finite wires if 
the screening from the substrate can be reduced. 

The central quantity of interest is the local density of 
states (LDOS) of inserting one particle at position x to 
reach an excited state {uj n | with Nq + 1 particles from the 
ground state |7V ) 

p(u,x) = ^2\(uj n \^ x \N )\ 2 S(uj-^ n ) 

n 

= Im / e iujt G R (t,x)dt. (1) 

n Jo 

Using the Luttinger liquid formalism it is now well un- 
derstood how the scaling laws of the LDOS and the cor- 
relation functions are affected by boundaries @, [Io|, E 
12, El Q Est The position resolved LDOS in finite 
wires was considered in Ref. [16| for the lowest few lev- 
els assuming a perfect degeneracy in the Luttinger liquid 
model. 

From the numerical side it is surprising that so far only 
very few works [13] were able to analyze the DOS in in- 
teracting lattice models exactly. Part of the difficulty 
appears to be that the LDOS in Eq. (pQ) is related to 
the time-dependent retarded Green's function G R (t). In 
dynamical DMRG methods dramatic progress has been 
made over the past few years. Recently, it was possi- 
ble to determine the spectral function as a function of 
the photoemission wavevector [17]. Unfortunately, indi- 



vidual levels have not yet been resolved due to a finite 
correlation time or a finite cutoff in the Lehmann rep- 
resentation of the Green's function. To our knowledge 
there have been no numerical calculations for the LDOS 
in interacting lattice models so far. 

In this paper, we are now able to quantitatively cal- 
culate the LDOS of a finite quantum wire by using a 
combination of bosonization and DMRG techniques [18[ . 
In the numerical DMRG calculations we follow a direct 
approach, by targeting several excited states (w n \ in the 
sector with one additional fermion on top of the ground 
state. Keeping track of the particle number and all the 
anticommuting transition operators it is then possible 
to evaluate the matrix elements of the LDOS in Eq. (pQ) 
with very good spatial and energy resolution compared 
to calculating the DOS via the Green's function. Using 
bosonization a general analytic formula for the LDOS of 
nearly degenerate states is derived. The combination of 
numerical and analytical calculations gives a good un- 
derstanding of the distribution of spectral weights over 
individual states. The well-known powerlaws, however, 
can only be observed after a summation over nearly de- 
generate states. The effect of spin is also discussed. 

We consider interacting spinless fermions hopping on 
a finite wire with L sites and "open" (^o = V'l+i = 0) 
boundary conditions 
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H = -tJ2 + ft-c.) + U n x n x +!, (2) 

x—1 x=l 

where n x =1^^:. This is the simplest model that il- 
lustrates Luttinger liquid physics for — 2t < U < 2t and 
is often used to describe quantum wires. It is believed 
that the results can be generalized to systems with spin 
as will be discussed later. Tunneling one single parti- 
cle (or hole) into the wire will create excited states that 
may in general involve several more fermions (i.e. a so- 
called "single particle" excitation may still be a many- 
body entangled state). In order to illustrate the na- 
ture of such many-body excitations, consider the system 
without interaction (U = 0) first. In that case, eigen- 
states are created from the vacuum by creation operators 
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FIG. 1: (color online) Single particle excitations on the rilled 
Fermi sea for the non-interacting system. 



lTT ^2 x sin(/c n + kp)x, where the wave vector 
k n = n-^j is measured relative to the Fermi point kp. 
For convenience kp = (No + is chosen so that 

the number of particles in the filled Fermi sea is given 
by No with energy Eq and there is no charging energy 
to the No + 1 particle sector as indicated in Fig. [TJ A 
typical example of a single-particle excitation in a tun- 
neling experiment is the creation of one fermion above 
the filled Fermi sea \a) = <4|iVo) as shown in Fig. [U 
In that case, the LDOS in Eq. ([T]) is simply given by 
the square of the corresponding standing wave. A sec- 
ond example of a possible excitation \b) = c\cqC-i\Nq) 
is shown in Fig. [TJ where also one additional fermion is 
created and another one has been excited from the Fermi 
sea. For non-interacting fermions the overlap matrix el- 
ements (b\ipl\No) of this many-body excited state vanish 
in the LDOS Eq. (pp), but this is no longer true for the in- 
teracting case. Using numerical results and bosonization 
we are now able to describe quantitatively how spectral 
weight is shifted from "simple excitations" (type \a)) to 
many-body excited states (type \b)) due to interactions 
and how this is reflected in the standing waves of the 
LDOS. 

Note, that if the dispersion relation was exactly linear, 
the states \a) and \b) would be degenerate. In fact a linear 
spectrum is often assumed, because in that case all many- 
body excitations are exactly quantized with energy levels 
— vpk n = j^Ti relative to Eq, where the number 
of possible partitions of n = Y^=i integers m t 

gives the degeneracy of the level uo n (e.g. for n = 3 = 2 + 
1 = 1 + 1 + 1 there are three possible partitions). However, 
numerically we keep track of the exact spectrum of the 
finite lattice model, which is never exactly linear so that 
the degeneracies are lifted. 

Even without assuming a linear spectrum the prob- 
lem can be bosonized by defining shifting operators g n = 
^£ c l+ n c £ f° r n 7^ 0? which obey bosonic commutation 
relations [g- n -,Qri] = nfin,n f on the infinitely continued 
spectrum [19]. The zero mode go = J2i ( c l c ^ — ( c £ c e)k F ^ 
counts the number of particles relative to the filled Fermi 
sea. It is possible to create any fermion state with a 



given particle number in terms of the g n . For example, 
we find by combinatorial methods that the addition of 
one fermion is represented by the following superposition 
of many-body bosonic states 



a*>- e ni(f)' 



|o>, 



(3) 



where |0) represents the ground state with Nq + 1 
fermions. Here, the sum runs over occupation numbers 
m i — (Q£Q-£)/£ which correspond to all possible parti- 
tions of n = ^2 £ ^m£ [20]. In particular, the states in 
Fig. [T] are given by 

\a/b) = 1 1 (q 2 1 ±q 2 )\0). (4) 

The local addition of one fermion is described by ipl ~ 

e~ ikFX ^ R (x) - e ikFX ^ R (-x), where 



^(x) =c(x)exp(^e- ife "^)exp(^ e - ifc "^) (5) 
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ik n x @r> 

n 



in accordance with open boundary conditions 0, [To|, E 
EH. Here, the prefactor c(x) = -j^—e^ ' 60 ^ also 

contains the zero mode operators £>o and 0o, which ex- 
actly create the (No + l)-particle ground state |0) = 
e i>(4>o-Qo z,+i)|jVo). The ambitious reader may enjoy ver- 
ifying that the state in Eq. Q is normalized and that 
the free fermion result | (A/o|c n ?/4| No) | 2 = j^-\sm(kp + 
k n )x\ 2 can also be evaluated in the boson language using 
Eq. © and ©. 

Interactions U can now be treated by a Bogoliubov 
transformation g n = ag n + (3g- n where the rescaling pa- 
rameters a = ±(1/VK + VK) and (3 = \(l/VK - Jk) 
are characterized by the Luttinger liquid parameter K. 
For the model in Eq. ([2]), K and the Fermi velocity 
vp are known exactly as a function of U by compar- 
ison with Bethe ansatz, e.g. at half filling 1/2K = 
1 - arccos(/7/2f)/7r and v F = 2tK sin(7r/2K)/(2K - 1). 
The interacting Hamiltonian becomes diagonal H = 
T+i Y^TLi Q-QtQ-t U P to non-linear corrections, which are 
described by higher order operators. After normal order- 
ing the prefactor of the transformed fermion operator in 
Eq. ([5]) is given by 



c(x) = 



cx(3 



(6) 



where the inconsequential zero modes are left out. Here 
a is a so-called cutoff parameter of order one lattice con- 
stant, which comes from a sum over physical states and 
can be determined numerically. 

The transformed boson vacuum and all eigenstates 
correspond to complicated many-fermion superpositions, 
which we can now analyze numerically. We implement 
the model in Eq. ^ with L = 78 sites and iV + 1 = 40 
using a multitarget DMRG. By keeping reflection and 
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FIG. 2: (Color online) DMRG data for the LDOS 
\(b\^(x)\N )\ 2 (first row) and \(a\^ (x)\N ) | 2 for L = 78, 
where \b) and |a) are eigenstates that have evolved from the 
n = 2 states in Fig. [T] when interactions C/ are switched on. 
The spectral weight is shifted towards \b) with increasing in- 
teractions, but | a) remains dominant. 



particle number symmetries, the position resolved ma- 
trix elements (E\^l,\No) of 67 excited states (E\ are cal- 
culated. The eigenenergies E agree to at least 8 digits 
with the exact Bethe ansatz. 

We use the simplest example of the states \a) and \b) 
in order to illustrate how the LDOS changes with U as 
shown in Fig. [2j The spectral weight of the many-body 
state | b) increases with interactions. The wavefunction 
can be analyzed with the help of Eq. ([5]), so that the 
eigenstates are known in terms of bosonic excitations, 
e.g. 1 6) « (0.541^ _ 0.455£ 2 )|0) for U = 1.6t. 

The exact bosonic states are determined by higher or- 
der operators in the bosonized Hamiltonian, which are 
not known quantitatively. However, since the complete 
nearly degenerate subspace for a given energy uo n = 
vpkn is known, it is in principle possible to calculate 
the sum of the LDOS of nearly degenerate states, i.e. 
Pn(x) = fupau p(w,x)du;. By summing over bosonic 
eigenstates Y\i=i (£ rn£r me\)~ 1 ^ 2 Q™* |0), corresponding to 
all partitions of n = X^ 777 ^ we nave found a general 
analytic expression for this "total LDOS" 



p n (x) = \c(x)\ 2 



E 



n 



c.c. 



(7) 



where Xi( x ) — OLe lkgX + (3e~ lke - x . This formula is the 
analytic generalization of the LDOS that previously was 
only known for the first three levels [16]. The total LDOS 
shows a pattern of standing waves that are modified by 
collective bosonic excitations in the form of characteristic 
density modulations as shown in Fig. [3l The analytical 
result in Eq. ([7|) can be compared with the numerical 
data as shown in Fig. [3] for the example of levels up to 
n = 5 for U = I At (K « 0.67). The overall scale is the 
only adjustable parameter, which determines the cutoff in 
Eq. ©, that ranges from a = 0.3 ... 0.7 as U = 0.2 ... l.St 




FIG. 3: (Color online) DMRG data for the total LDOS of 
different levels n for U = 1.4t (blue crosses) compared to the 
bosonization results in Eq. © with a = 0.56 (solid line). Red 
dots mark the LDOS from bosonization at discrete lattice 
points. 



is increased. Note, that the lattice structure of the model 
([2]) is incommensurate with the rapid oscillations of the 
LDOS, which leads to a beating of the signal. The quan- 
titative agreement with the analytical prediction is very 
good up to small deviations near the boundary. 

The number of nearly degenerate states at each en- 
ergy uo n grows rapidly with the number of partitions 
~ exp(7r ^2n/?>) /4n>/3 and each of those states has an 
increasingly complex wavefunction, so that a description 
in terms of approximate power laws is normally used. In- 
deed, if the position dependence is integrated over, the 
total DOS p n = J2 X Pn( x ) m Eq. ([7j) should follow an 
approximate powerlaw 



1 / 7ra 

T(l+2/3 2 ) y L+ 
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r(i+2/3 2 ) 
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2(3 2 



(8) 



However, this only holds for the sum over all nearly de- 
generate states. The individual states on the other hand 
show a much more interesting energy dependence of the 
position integrated DOS as shown in Fig. [H The domi- 
nant spectral weight in each multiplet has evolved from 
the original simple excitation in Eq. (|3]) (type \a)), which 
has a surprisingly flat energy dependence. The many- 
fermion states (type |6)) are much more numerous, but 
have a much smaller DOS. 

The position integrated total DOS of each level in Fig. [4] 
agrees well with the integral of Eq. ([7j) . There is no sign of 
any deviations from the low energy theory, which should 
break down as the cutoff is approached uo ~ vp/a, i.e. be- 
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FIG. 4: (Color online) DMRG data for U = 1M and L = 78 
of the position integrated DOS p(oS) — p(uj,x) for each 
individual state (blue crosses). The green stars represent the 
sum over all nearly degenerate states at each discrete energy 
cj n , which agrees very well with the integral of Eq. J7J) (cir- 
cles) . 



fore the DOS reaches the non-interacting value of unity. 
The alternation of the analytical results relative to the 
powerlaw is due to boundary effects, which induce addi- 
tional oscillations with sh\2ux/ vp [14]]. In general, for 
longer range interactions it is also possible to use a mo- 
mentum dependent Luttinger liquid parameter [21]. 

The situation as shown in Fig. [4] is actually quite 
generic for interacting systems also in higher dimensions: 
The single particle spectral weight is spread over many- 
body excitations, thereby renormalizing position, width, 
and total weight. Our results show quantitatively how 
an effective finite lifetime of quasiparticle states evolves 
due to a spread over many-body states. 

In one-dimensional systems with spin degrees of free- 
dom (7 =|, j, the excitations can be described by two 
independent bosonic theories for spin and charge [191 ]. 
The spectrum is again given in the form of nearly de- 
generate multiplets uo = uo ns + uo Uc , but now labeled by 
two quantum numbers (n s , n c ) for spin and charge with 
^n s = v s k ns and co nc = v c k nc1 where v c = v s in the non- 
interacting case. A single particle state a \No) is now 
divided over all spin and charge multiplets (n s , n c ) cor- 
responding to the different possibilities to write the sum 
n = n s + n c . However, for a given n = n s + n c the multi- 
plets (n s , n c ) may now be of quite different energy, since 
in general v c > v s for repulsive interactions. The spec- 
tral weight is therefore spread into several non- degenerate 
spin and charge multiplets, which is fundamentally dif- 
ferent from higher dimensional systems. Moreover, the 
degeneracy within each multiplet is also lifted in the same 
fashion as above. In each multiplet there is again exactly 
one dominant state. The LDOS also shows a very in- 



teresting superposition of spin and charge density waves 
[l6[ , which can be calculated in detail with rescaled waves 
Xi{x)/y/2 in Eq. ©. However, the smoking gun for Lut- 
tinger liquid behavior in a finite wire would be to identify 
the increasing number of states due to spin-charge sep- 
aration and the degeneracy splitting with increasing a;, 
which is very different from an equally spaced single par- 
ticle spectrum. Numerically, however, the important ef- 
fects of interactions can be much better illustrated using 
the single channel (spinless) case discussed here, which 
avoids the complications of a huge number of states and 
a superposition of spin and charge standing waves. 

In summary, we were able to describe the LDOS for fi- 
nite quantum wires in detail by a combination of analyt- 
ical and numerical methods. The analytical results allow 
the calculation of the total LDOS in each multiplet. The 
DMRG calculations give the detailed distribution of lo- 
cal spectral weights over all many-body states in the low 
energy region. The combination of both methods shows 
explicitly how the wavefunctions and boson representa- 
tions of excited states evolve as interactions are turned 
on. The cutoff parameter a(U) and the normalization in 
Eq. (|6]) can be determined from DMRG for the standard 
model in Eq. ([2]). Our results show that powerlaws are 
not sufficient to adequately describe the low energy be- 
havior of individual levels. Instead, a large number of 
discrete states with varying spectral weights and oscil- 
lating wavefunctions would be the generic signature of 
Luttinger liquid behavior in finite wires. The number of 
many-body states with non-zero spectral weight is small 
for the first few levels, but increases exponentially with 
energy. 
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